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Risk Extrapolation for Chlorinated Methanes as Promoters vs Initiators of Multistage Carci¬ 
nogenesis. Bogen. K. T. ( 1990). Fnndam. Appl. Toyicat. 15, 536-557. “Cell-kinetic multi¬ 
stage" (CKM) models account for clonal growth of intermediate, premaiignant cel! populations 
and thus distinguish somatic mutations and cell proliferation as separate processes that may 
influence observed rates o! tumor formation. This paper illustrates the application of two ver¬ 
sions of a two-stage CKM model 1 one assuming exponential and the other geometric prolifera¬ 
tion ol intermediate ceils) for extrapolating cancer risk potentially associated with exposure 
to carbon tetrachloride, chloroform, and dichloromethane, three suspect human carcinogens 
commonly present in trace amounts in groundwater supplies used for domestic consumption. 
For each compound, the models were used to calculate a daily orai “virtually safe dose" (VSD) 
to humans associated with a cancer risk of I0' 5 . extrapolated from bioassay data on increased 
hepatocellular tumor incidence in B6C3FI mice. Exposure-induced bioassay tumor responses 
were assumed iirst to be due solely to "promotion" (enhanced proliferation of premaiignant 
cells, here associated with cytotoxicity), in accordance with the majority of available data on in 
vivo genotoxicity for these compounds. Available data were used to model dose response for 
induced hepatocellular proliferation in mice for each compound. Physiologically based pharma¬ 
cokinetic models were used to predict the hepatotoxic effective [ metabolized) dose as a function 
of parent compound administered dose tn mice and m humans. Resulting calculated VSDs are 
shown to be from three to five orders of magnitude greater than corresponding values obtained 
assuming each of the compounds is carcinogenic only through induced somatic mutations 
within the CKM framework. Key issues and uncertainties in applying CKM models to risk 
assessment for cancer promoters are discussed, o i wo Society of Toxicology. 
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Multistage cancer risk models that account 
for clonal growth of intermediate (initiated, 
premaiignant) cell populations, or “cell-ki¬ 
netic multistage” (CKM) models, have long 
been successfully used to describe age-specific 
incidence trends for human cancer and ex¬ 
perimentally induced tumors in animals 
(Armitage and Doll, 1957, 1961). More re- 


1 To whom ail Correspondence should be addressed at 
Environmental Sciences Division (L-453). Lawrence 
Livermore National Laboratory. University of Califor¬ 
nia. P.O. Box 80S. Livermore. CA 94550. 


cently, a two-stage stochastic model allowing 
for birth and death of initial and intermediate 
cells has been used to mode! human cancer 
incidence and has also been proposed to de¬ 
scribe spontaneous and induced carcinogene¬ 
sis in general (Moolgavkar and Venzon, 
1979; Moolgavkar et a!,, 1980; Moolgavkar 
and Knudson, 1981). There is growing inter¬ 
est in applying such a model to the problem 
of predicting human cancer risk posed by ex¬ 
posure to environmental agents (Mool¬ 
gavkar, 1983; Thorslund et al., 1987; Swen- 
berg et ai, 1987; Moolgavkar et ai., 1988; 
Wilson, 1988; Dewanji et al, 1989: Bogen, 
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1989). Because these models distinguish be¬ 
tween mutations and cell population kinetics 
as separate processes that influence the rate 
of tumor formation, such models can predict 
widely different cancer risks for agents affect¬ 
ing only one process (e.g., “promotion” via 
induced net cell proliferation assum¬ 
ing a threshold- or quasithreshold-type dose 
response) but not the other(e.g., somatic mu- 
tatioin assuming a no-threshold dose re¬ 
sponse). Such an approach allows for distinct 
treatment of genotoxic carcinogens and non- 
genotoxic tumor promoters, provided that 
generally acceptable empirical criteria For the 
latter classification can be specified. This ap¬ 
proach has been advocated (Williams. 1967; 
Weisburger and Williams, 1983; Wilson, 
1988) but is currently not specifically incor¬ 
porated into widely used procedures for can¬ 
cer risk assessment (U.S. EPA, 19861. 

Specific proposals for incorporating infor¬ 
mation on induced cell proliferation into 
quantitative risk assessment for cancer pro¬ 
moters currently focus on the use of tradi¬ 
tionally defined no-observed-effect levels and 
safety factors (Weisburger, 1988; Clayson et 
at.. 1989; Butterworth, 1989). In contrast, 
this paper illustrates a more rigorous ap¬ 
proach whereby an approximate two-stage 
CKM model 5 (Armitage and Doll, 1957; 
Moolgavkar, 1983; Thorslund etal. 1987) is 
incorporated along with experimental data 
on tumor induction and induced cellular pro- 

2 The approximate CKM model referred to here is 
based on a “semistochastic” two-stage CKM model in 
which interstage transition via somatic m utation is mod¬ 
eled as a stochastic process, but in which proliferation 
of premalignant cel] clones is modeled as deterministic 
exponential growth. In the “fully stochastic” two-stage 
CKM model proposed by Moolgavkar et a!. (see refer¬ 
ences cited), both the transition and proliferation pro¬ 
cesses are modeled as stochastic birth-and-death pro¬ 
cesses. Approximate forms of the semistoehastic model 
are considered here to illustrate the application of CKM 
models to quantitative risk assessment for cancer pro¬ 
moters. However, these approximations are identical to 
analogous approximations for the corresponding “fully 
stochastic" model (see Appendix.!. 


liferalion into a risk extrapolation procedure 
for promoters of hepatocellular carcinoma in 
mice. Because this approximate CKM model 
assumes that premalignant cell clones grow 
exponentially in vivo, a new version of this 
model (Bogen. 1989), using the biologically 
more plausible assumption of geometric cell 
proliferation, was also applied in the present 
analysis. The risk extrapolation procedure 
was applied to the chlorinated methanes, car¬ 
bon tetrachloride (CCI 4 ), chloroform (CHCT), 
and dichloromethane (methylene chloride, 
CH;Clj), three compounds known to be hep- 
atocarcinogenic to mice and also known (as 
discussed below) to be inefficient or ineffec¬ 
tive at inducing genotoxicity in mammalian 
cells in vivo. 

In the following analysis, all three com¬ 
pounds are first assumed to be “pure promot¬ 
ers” within theCKM framework, an assump¬ 
tion that is suggested (but by no means 
proven) by data indicating minimal or non¬ 
existent genotoxicity of these compounds in 
mammalian cells in vivo coupled with an 
ability to induce hepatocellular toxicity and 
proliferation (discussed below). Thus, in the 
first part of this analysis, the compounds are 
presumed ineffective at producing mutations 
assumed to be necessary in a CKM carcino¬ 
genic process; rather, they are assumed to be 
effective, at cytotoxic doses, in selectively in¬ 
creasing the rate of net cellular proliferation 
within premalignant cell clones. Available 
data are used to model the dose-response re¬ 
lationship for compound-induced cell prolif¬ 
eration in mice, and this relationship is as¬ 
sumed to be the same for humans. Increased 
predicted cell proliferation is then presumed 
to be proportional to increased growth pa¬ 
rameter values in the CKM models fit to cor¬ 
responding sets of mouse bioassay data on 
liver tumors induced by these compounds. For 
each compound, CKM models are then used to 
extrapolate a corresponding “virtually safe 
dose” (VSD) associated with a predicted in¬ 
creased cancer risk to humans of 10 -6 . 
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In the second part of this analysis, the 
VSDs calculated using CKM “promotion” 
models are compared to values calculated us¬ 
ing corresponding CKM models under the 
assumption that the compounds considered 
are purely mutagens (and not promoters) in 
the CKM context. In applying both the CKM 
promotion-only and mutation-only models, 
appropriate physiologically based pharmaco¬ 
kinetic (PBPK) models are used throughout 
to predict the cytotoxic or genotoxic effective 
(metabolized> dose to liver in mice and in hu¬ 
mans as a function of applied dose of parent 
compound. Finally, the calculated CKM- 
based VSDsare compared to VSDs based on 
the same mouse bioassay data sets but calcu¬ 
lated using the “linearized” multistage model 
currently used in many regulatory cancer risk 
assessments, and issues concerning the appli¬ 
cation of the outlined procedure to cancer 
risk assessment are discussed. 


Hepatocarcinogeniciiy of Chlorinated Meth¬ 
anes in B6C3F1 Mice 

Bioassay data. Experimentally induced 
caicinogemc effects associated with lifetime 
exposure to CCU, CHC1 3 , and CH 2 CI 2 in ro¬ 
dents have implicated these compounds as 
probable or potential human carcinogens 
(IARC, 1979a,b, 1982, 1986; U.S. EPA, 
1984, I985a,b,c, 1987). The present CKM- 
based cancer risk extrapolation was based on 
bioassay data on hepatocarcinomas (plus ad¬ 
enomas for CH 2 CF) in control and low-dose 
female B6C3F1 mice in studies by the Na¬ 
tional Cancer Institute (NCI, 1976a,b, 1977; 
Reuber, 1979) and the National Toxicology 
Program (NTP, 1986), as summarized in Ta¬ 
ble 1. With regard to the bioassay data for 
CC1 4 , note that because 100% incidence was 
observed in the 40 low-dose female mice used 
and thus lower dose rates might also have 
yielded a 100% incidence, these data cannot 
be used to derive a point estimate of hepato- 
carcinogemc potency for CCU in mice. Nev¬ 


ertheless, the data suggest that an oral dose of 
1250 mg/kg CC1 4 on the NCI dosing regimen 
would yield a cumulative liver cancer risk 
greater than (40- !)/40or97.5%by97 weeks 
of age in these mice. Under this assumption, 
lower-bound CKM-based risk extrapolations 
were undertaken for CCU based on the NCI 
mouse bioassay data. 

Information Concerning Likely Mechanism 

of Carcinogenic Action 

CCI 4 . Most assays of CCl 4 ’s mutagenicity 
have been negative, particularly those in vol v- j 
ing mammalian systems (U.S. EPA. 1984), | 

In vivo assays of hepatic DNA damage using 7 
rats exposed to CC1 4 have been uniformly 
negative (Mirsaiis and Butterworth, 1980; j 
Mirsalis et a/.. 1982; Bermudez eta/.. 1982). , 

In contrast, CCU is a classic agent for indue- 1 
ing hepatic necrosis and regenerative hepato¬ 
cellular proliferation (Leevyef al.. 1959; Din- j 
man et ai, 1963; Scbultze et al.. 1973) and ' 
has been shown to induce in exposed mice a 1 
substantia] increase in the percentage of he- 
patocyles in S-phase compared to that in un¬ 
exposed mice (Mirsalis et al.. 1985; Doolittle 
et al.. 1987). Chemically induced hepatocar- 
cinogenesis has been shown to be promoted 
by necrogenic doses of CCU in B6C3F1 mice j 
and in hamsters (Dragani et al. 1986; Ta¬ 
naka et ai, 1987); in mice this has been I 
shown to be associated with increased general j 
hepatocellular proliferation plus decreased 
cytotoxicity specifically to neoplastic hepato- j 

cellular nodules, i.e., with increased net pro- | 
liferation of premalignant hepatocyte clones 
(Dragani et al.. 1986). A number of effective 
tumor-promoting compounds, including 1 

phorbol esters and CC1 4 , have also been 
shown to be potent activators of intracellular 
protein kinase C, an enzyme thought to be 
involved in regulating cell growth (Rogham 
etal., 1987). 

CHClj . Most genotoxicity assays for 
CHClj have yielded either negative or incon- 
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TABLEI 

Input Parameter Values for CKM Models opHepatocellulah Carcinoma 0 Incidence 
in Mice Exposed to Chlorinated Methanes 




Estimated 
effective dose, d c 
(mg/kg liver) 



Input parameter values 

«/ 


Compound 

Applied dose'’ 

t 1 (year) 

(2 (year) 

/, (year) 

R„(%) 

RA%) 

Xp(year 2 
X tO' 6 ) 

CCI 4 

1250 

21.1 

7 

52 

85 

52 

97 

52 

1,25 

>97.5 

1500 

CHC1, 

230 

2230. 

52 

& 

52 

s 

1.25' 

88.9 

1500 

CH.Cl, 

2000 

1210 

hi 

52 

im 

51 

112.5 

52 

6.0 

34.8 

1750 


a Hepatocellular carcinoma or adenoma in the case ofCH.CL. 

* CC1 4 , dose I in mg/kg) to 28-g female B6C3FI mice was administered orally 5 days/week for 78 weeks in com oil 
vehicle (NCI, 1976a, 1976b. 1977). CHCl). dose (in mg/kg) to 30-g female B6C3F1 mice was administered orally 5 
days/week for 78 weeks in com oil vehicle (NCI, 1976a: Reuber, 1979). CH 2 C1,. dose (in ppm) to 32-g female B6C3FI 
mice was administered in respired air 6 hr/day, 5 days/week for 102 weeks 1NTP. 1986). 

'Toxicology effective dose to hepatocytes (total metabolized dose 1 hr postexposure for CC1 4 and CHC1 3 , total 
daily GST-metabolized dose to liver for CH 2 C1 3 ) was estimated using physiologically based pharmacokinetic models 
described in the text. 

d See text for explanation of CKM dose-response models and their input parameters. 

’'Control hepatocellular carcinoma incidence was0/19: R 0 was assumed to be equal to the value pertaining to CCt 4 . 


elusive results (Bridges et at., 1981; Brookes 
and Preston. 1981; U.S. EPA 1985a; Rosen¬ 
thal, 1987). In contrast, hepatotoxicity has 
been observed to occur in male and female 
B6C3F1 mice chronically exposed to oral 
doses of CHCI 3 in corn oil or in water (Bull 
et a/., 1986). Hepatocellular regeneration 
was found to be induced by this compound 
in liver and kidney of male B6C3F1 mice ex¬ 
posed to single oral doses, as measured by rel¬ 
ative DNA incorporation of [ 3 H]thymidine 
(Reitz etal., 1982). A similar finding was ob¬ 
served for male Swiss albino mice exposed to 
single oral doses of chloroform in com oil and 
toothpaste vehicles { Moore etal., 1982). Per¬ 
eira et at. (1984) also found replicative hepa¬ 
tocellular DNA synthesis, measured by 
[ 3 H]thymidine incorporation into DNA, to 
be enhanced in male B6C3F1 mice after ex¬ 
posure to single oral CHC1 3 doses. Like CCf, 
CHC1 3 has been shown to be a potent activa¬ 
tor of protein kinase C (Roghani et a!., 
1987). 


Cff 2 Cl z . With the exception of a few re¬ 
pons of induced chromosome damage in vi¬ 
tro (Jongen et aL, 1981; Thilagar and Ku- 
maroo. 1983: Thilagar et at.. 1984). this com¬ 
pound has generally not been found to be 
genotoxic to mammalian cells in vitro 
(Gocke el a!., 1981; Perocco and Prodi, 1981; 
Jongen et aL. 1981; Andrae and Wolff, 1983; 
Burek el aL, 1984). It also has been repeat¬ 
edly shown to be negative in in vivo assays of 
DNA damage (Burek et at ., 1984; Sheldon et 
aL, 1987: Trueman and Ashby, 1987; Green 
el aL, 1988; Raje et at., 1988; Ottenwalder 
and Peter. 1989), except for one recent report 
of slightly increased DNA damage measured 
by alkaline elution of DNA from hepatocytes 
of rats orally exposed to 2.55 g/kg CH 2 C1 2 in 
com oil over 24 hr (Kitchin and Brown, 
1989). In contrast, chronic exposure of rats 
or mice to this compound at levels as low as 
100 to 500 ppm by inhalation or 170 (rats) 
to 590 (mice) mg/kg/day by ingestion in 
drinking water is known to cause mild hepa- 
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tocellular and renal toxicity (Klassen and 
Plaa, 1966; Weinstein etal., 1972; U.S. EPA, 
1985b; Kirschman et al. 1986; Nitschke et 
al, 1988), although CHjCL is known to be 
much less hepatotoxic than CC4, CHC! 3 , 
and several chlorinated ethanes (Klassen and 
Plaa, 1966). Although chronic exposure of 
female Sprague-Dawley rats to concentra¬ 
tions up to 500 ppm in air resulted in no evi¬ 
dence of increased hepatocellular DNA syn¬ 
thesis (Nitschke et al., 1988), single inhala¬ 
tion exposures of B6C3F1 mice to 4000 ppm 
CH 2 C4 did result in a small but statistically 
significant increase in the percentage of hepa- 
tocytes in S-pbase in these animals (LeFevre 
and Ashby, 1989). Like CCl 4 and CHClj, 
CH 2 CL has been shown to be a potent activa¬ 
tor of protein kinase C fRogham el al., 
1987). 


METHODS 


PBPK Models 


CC1 4 . Metabolism of this compound is known to be 
the cause of its hepatotoxicity in rats (Reynolds. 1967; 
Bruckner et al- 1986: Shertzer el al.. 1987. l988;Labbe 
! 't al. 1987) a-d ir, one study the rate of production of 
CO, (the primarily metabolite in rats) within 1 hr of dos¬ 
ing was shown to be the measure of CC1 4 metabolism 
most highly correlated ( R 2 = 0.88) with the observed de¬ 
gree of enzyme-assayed liver injury (Reynolds el al.. 
1984). Therefore, in the present analysis, the total 
amount of CC1 4 metabolized in 1 hr posiexposure was 
assumed to be an appropriate measure of hepatotoxic 
effective dose. To estimate uptake and metabolism of 
CCU in mice and humans, the PBPK model of Pausten¬ 
bach el al. ( 1988), based on the four-compartment 
model developed by Ramsey and Andersen (1984). was 
adapted with the following modifications. CCI 4 metabo- 
[ism in rats was assumed to occur by a single saturable 
palhway with equal to 0.14 mg/hr for a 225-g rat 
(Gargas el al.. 1986a). (Corresponding F mlI values 
scaled, usmg the methodology of Paustenbach el al. 

(1988), for modeling 30-g mice and 70-kg humans were 
0.034 and 7.78 mg/hr, respectively.) Blood/air and 
tissue/air partition coefficients for rats and humans used 
were those of Paustenbach el al. ( 1988); those for mice 
were set equal to those for rats. Alt other physiological 
parameter values were taken to be the reference values 
reported by the U.S. EPA (1988). except that the latter 


reported alveolar ventilation rate for 70-kg humans was 
increased from 5.0 to a more realistic 5.6 liters/min. 
PBPK parameter values for mice and humans with 
weights different from those above were scaled by the 
method described by Paustenbach et al. ( 1988). 

A study by Withey ei at- ( 1983) of the absorption ki¬ 
netics of four chlorinated organic compounds adminis¬ 
tered by a single gavagedosc to rats in aqueous and com 
oil vehicles demonstrated that, over a 5-hr observation 
period, a corn oil vehicle decreases the rate and extent of 
uptake of these compounds relative lo that observed us¬ 
ing an aqueous vehicle; the greatest decrease in uptake 
observed was for trichloroethylene, the least water-solu¬ 
ble compound tested. In the present analysis, oral ab¬ 
sorption kinetics for CC1 4 were presumed to be similar to 
those observed by Withey etal (1983) for trichloroethyl¬ 
ene administered in 10 ml/kg of com oil to 400-g Wistar 
rats, which yielded a roughly constant, 90-mtn absorp¬ 
tive pulse equivalent to a clearance of compound from 
the vehicle at a rate of about 0.11 ml vehicle cleared/kg 
body wt/min, resulting in a total area under rhe blood 
concentration umestimecurvetor AUBC. here assumed 
to be proportional to total absorbed dose) of about one 
170th of that obtained using an aqueous vehicle (here 
assumed to reflect 100% uptake of the administered 
dose). This vehicle clearance or ingestivc uptake rate 
(times 1 70~ 1 ) was assumed to Scale to the -0.36 power 
of body weight 1 Bogen, 1988) and to apply only to single¬ 
dose (and not to chronic) exposure scenarios involving 
a corn oil vehicle. The effect of such a vehicle on absorp¬ 
tion efficiency of chronically administered oral doses of 
CCI 4 and similar compounds is currently unknown. 

CHClj. Chloroform is metabolized m the tivei to a re¬ 
active intermediate (phosgene) which is rhe proximate 
cause of this compound’s hepatotoxicity (llett el al. 
1973; Brown et al.. 1974a; Lavigne and Marchand, 1974; 
McMartin ei al, 1981; U.S. EPA, 1985a). To estimate 
effective (metabolized) dose as a function of applied 
dose, a model of CHClj metabolism was used based on 
the following CHCl^-specific parameterization of the 
PBPK model of Ramsey and Andersen (1984). Tissue/ 
air and blood/air partition coefficients used for mice 
were assumed to be equal to those obtained for rats and 
that obtained tor mice, respectively, by Gargas et al 
(198S. 1989) by in viiro techniques. Physiological pa¬ 
rameter values used for mice ant) humans were the refer¬ 
ence values cited above (U.S. EPA, 1988). except that 
the alveolar ventilation rate for 70-kg humans was in¬ 
creased to 378 liters/hr, based on the data reported in 
U.S. EPA (1988). The blood/airpariilion coefficient for 
humans was taken to be 8.0. the median of several re¬ 
ported values(Steward et at .. 1973; Sato and Nakajima. 
1979; Gargas et al. 1989). Human tissue/blood parti¬ 
tion coefficients were similarly taken from the literature; 
2.1 for liver and 35 for fat (Steward et al. 1973). 1.9 for 
the richly perfused compartment (Fcingold and Hola- 
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day, 1977). and 1.5 for the muscle/skin compartment 
(average of values for “muscle” and “vessel-poor'' com¬ 
partments reported by Feirtgold and Holaday, 1977). 
The Michaelis constant for mice and humans was set 
equal to0.25 mg/liter, a value measured for ratsfCargas 
eiai., 1988). 

Uptake of single doses of ingested CHC1 3 in com oil by 
rats was observed by Withev el ai (1983) to yield a 
roughly constant pulse of absorption into blood lasting 
about 30 min. This corresponds to a rate of uptake from 
vehicle of 0.33 ml vehicle cleared/kg body wt/min. 
which was assumed in the present analysis to scale to the 
—0.26 power of body weight and to apply only to single¬ 
dose exposures involving a corn oil vehicle. Mass-bal¬ 
ance metabolism studies ofCHCIi orally administered in 
corn oil to mice and rats have shown that the majority of 
such a dose is metabolized to CO, within a period of 
hours after administration (Reynolds et ai. 1984; Mink 
ei at.. 1986). Accordingly, it was assumed that mice ab¬ 
sorb lOOTf of orally applied doses of CCh in corn oil. and 
that, as with CC1 4 . the total amount of CHClj metabo¬ 
lized in 1 hr postexposure is the appropriate measure of 
the hepatotoxic effective dose. 

f ma >. for 25-g mice (4.7 mg/hr) was estimated by fit¬ 
ting the metabolized dose predicted by the PBPK model, 
parameterized as stated except for T„„.to the average of 
the metabolized amounts observed for mice given a sin¬ 
gle gavage dose of 60 nig/kg CClj in corn oil ( Brown el 
ai. 1974b; Taylor et ai.. 1974). l' m „ for a reference 70- 
kg human (560 mg/hr) was similarly estimated using 
data on metabolism, collected over 8 hr. from human 
volunteers given single oral doses of 500 mg CClj in olive 
Oil (Fry t! ai.. 1572). FEiPK. paittntetci values for mice 
and humans with different weights were scaled as de¬ 
scribed above. 

CHiClz- Metabolism of this compound yields reactive 
species through both a saturable (oxidaiive 
mitochondrial) mixed-function oxidase pathway and a 
relatively nonsaturable (nonoxidative cytosolic) gluta¬ 
thione 5-transferase (GST) pathway (Kubic el ai.. 1974: 
Ahmed el ai.. 1980; Gargas el at.. 1986b; Reitz el ai.. 
1986). These reactive metabolites, as those from other 
multihalogenated methanes, are likely to be responsible 
for exposure-related hepatotoxieity because the parent 
compound is relatively nonreactive with tissue macro- 
molecules i Ahmed et ai., 1980). Consistent with this hy¬ 
pothesis is the similarity between hepatic eytoxicity ob¬ 
served in mice given acuteorsubehronicdosesofCHjCl, 
and that observed in mice treated with CCLj (Weinstein 
et at. 1972). In the present analysis, hepatotoxieity was 
presumed to arise only from that fraction of daily applied 
dose metabolized in liver by the GST pathway, which 
is consistent with the observed pattern of experimental 
CHjCJj- induced mutagenicity in bacteria and lumorigc- 
nicity in rodents (Green et ai., 1988; Andersen ei ai. 
1987; U.S. EPA. 1987). To estimate this effective dose 


asa function ofapplieddose.lheRritzeraf. (1988,1989) 
parameterization of the PBPK model of Andersen el ai 
(1987) for CH 2 C! 2 in mice and humans was used. PBPK 
parameter values for mice and humans with different 
weights were scaled as previously described. 


Mode! for induced Celt Proliferation 

Dose-response data on induced hepatocellular prolif¬ 
eration in mice for the three compounds considered were 
obtained from the studies listed in Table 2. To make use 
of the oral exposure studies involving CHClj and CCI 4 , 
oral uptake efficiency was estimated as described above. 
In the case ofCCI 4 . different groups ofthe30-g mice used 
in the Doolittle ei ai. (] 987) study were given single oral 
doses of from 0 to 100 mg/kg CCI 4 in 10 ml/kg corn oil 
for different study periods ranging from I to 14 days. For 
the present analysis, metabolized doses. ,0. correspond¬ 
ing to these applied daily doses were estimated as I / 

170th of the ingested CC1 4 , cleared from the vehicle at a 
constant rate of 0.22 ml/kg/min for 45 min, and in¬ 
jected into the liver compartment. Similarly, it was as¬ 
sumed for the present analysis that the oral doses of chlo¬ 
roform in the Reitz ei at ( 1982) study were administered 
to 30-gmicein lOml/kg com oil vehicle, such that 1 00 % 
uptake in the dosed mice occurred by constant infusion 
into the liver compartment over a 15-min period. The 
question of oral uptake efficiency does not arise in the 
ease of the LeFevre and Ashby (1989) study of hepato¬ 
cellular cell proliferation induced by CH 2 CF delivered by 
inhalation. 

The relationship between estimated values of acute 
effective dose D and corresponding predicted values of 
the percentage of hepalocytes in S-phase (?o5) was mod¬ 
eled as a threshold-like relationship based on:( !) the ob¬ 
vious nonlinearity of the corresponding dose-response 
data for CHCI, and CC1 4 , the two compounds for which 
sufficient relevant data were available to make such a di¬ 
rect determination (see Table 2, Fig. 1); (2) a specific 
assumption made here that such proliferation represents 
regenerative ceil growth following necrosis caused by 
chemically induced cytotoxicity; and (3) the traditional 
use of threshold models for virtually all endpoints other 
than cancer and mutation in bolh descriptive and regula¬ 
tory toxicology. Clearly, apparent nonlinearity in ob¬ 
served dose-response data can never disprove a linear 
dose-response relationship at very low doses for which 
data are not available. Furthermore, the choice of a 
threshold-like model is critical in the present analysts, be¬ 
cause the use of a nonthreshold relationship (such as one 
including a linear function Of effective dose) here would 
result in dramatically higher CKM-derived estimates of 
predicted increased risk over background associated with 
exposure to low doses (see Discussion). Therefore, the 
rationale for using a threshold-type dose response for in- 
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TABLE 2 


Lognormal Model of Hepatocellular Proliferation in Mice Exposed to Chlorinated Methanes 


Compound 

Data source: strain; 
sex: weight 

Assay type 0 

Applied 

dose* 

Estimated 
effective dose*". 

D{ mg/kg liver) 

%S (%) 

Estimated parameters of dose- 
response model d 

«(8/kg 

liver) a p 

o 

o 

Doolittle ft ill. 

^5", vtr 

0 <61 


0.5 



(IUS7): 


1 (fVl 

0 

3.0 



CD-I: M:30g 


2'(3) 

0.046 

1.0 





in (6i 

(1.092 

S.O 





20'H) 

0.45 

2.0 





25 (6) 

1.1 

2.0 

0.0017 1.25 1.00 




30 (6) 

1.3 

5.5 





40 16) 

1.7 

19 





50 161 

2 i 

22.5 

1 




50 (31 

j 

30 





100*' (3) 

3.9 

34 

1 




100 ft) 

3.9 

27 

- 

CHC'lj 

Reitz vl ill. ()9&2): 

‘ll-DNA. vi\ 

0 1C) 

0 

0,5* 

| 


BftOFl: M: 


15 <5) 

27 i 

0.5 

4.03 2.02 0.20 


30 g ; 


60 i‘6> 

937 

1.14 

t 




240 (6) 

2740 

7.23 


CH,Cl, 

L dev re and AsJ‘ih> 

' r$. V|\ 

t)(5) 

0 

0,026 

5.79 2.02* 0.013! 1 


if989); B6C3FI; 


4000/5) 

902 

0.162 

1 


M; 26 g 





! 


determination of percentage of hepatocytes in S-phase i' e ”-5) isolated 48 (CCI 4 ) or 24 (CH,CL) hr after oral 
administration of compound or vehicle alone; nr. in vitro ireatment of hepatocytes with [ 3 H)thymidine ( J H-T) after 
isolation by liver perlusion: viv, ip injection of animals with J H-T prior to isolation of hepatocytes by liver perfusion. 

H-DNa, cellular regeneration index estimated by determining relative incorporation of J H-T in liver DNA. See 
corresponding references for details. 

L-.... [p mg/kg oiaily aduiiiustered in com oil vehicle (CCL, CHC1 3 ) or in ppm administered in respired air 
(CH.CJ,): number of animals used is shown in parentheses. 

‘ Toxicologicallv effective dose to hepatocytes (total metabolised dose 5-hr postexposurc for CCL and CHCL, 
total GS (-metabolized dose to liver for CH;CL } was estimated using physiologically based pharmacokinetic models 
described in the text. 

d See text for explanation oflognormal dose-response model and its parameters. %Sn, %.S for D = 0. Values for a 
and <r are least-squares estimates based on the model fit to the corresponding data: values for » are based on and 
corresponding values of effective bioassay dosed from Table I. as explained in the text. 

'Three animals received the dose at TimeO and 24 hr (total oftwodoses). 

' Assumed weigh! (actual value unspecified). 

'Data not reported as %S. but rather as 3 H-T in DNA. %£„ assumed to be equal to the value pertaining to CC1.; 
other %S values were assumed to be equal to ffeSo times the ratio of observed levels of J H-T in DNA of dosed animals 
to that in controls. 

h Insufficient data for estimation ol both w and o'. o was assumed to be equal to the value pertaining to CHClj. 


duced hepatocellular proliferation will be discussed in 
detail. 

The modeling choice that was made here reflects the 
difference in probable causal mechanism between in¬ 
duced ceil proliferation and toxic endpoints that might 
reasonably be associated with a linear dose-response at 
low doses. The discrete, quantal, and potentially irrevers¬ 
ible nature of DNA mutation caused by exposure to cer¬ 
tain reactive chemicals or ionizing radiations is both the¬ 


oretically and experimentally associated with a linear 
nonthreshold dose-response relationship at very low 
doses (NAS. 1980; NCRP, 1989). Indeed, this fact has 
been the basis for the traditional application of linear 
nonthreshold models in risk prediction for environmen¬ 
tal carcinogens (U.S. EPA. 1980; Albert, 1981 )and is the 
basis for the present application of a linear function of 
effective dose in the mutation component of the CKM 
models considered here (discussed below). 
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rs of dose- 
iei* 



1 *hr 30sa (mg/kg livery 

Fig. 1. Lognormal model (see Table 2 and corre¬ 
sponding text) fn to the data of Doolittle et at. ( 1987) on 
the percentage of hepatocytes in S-phase (%S) in male 
CD-I mice 48 hr after exposure to different doses of CCU 
in corn oil. administered orally either once or twice sepa¬ 
rated by 24 hr (see Table 2). 


It is certainly conceivable that a single molecular inter¬ 
action might be capable of causing cell replication either 
directly, or indirectly by causing cell death triggering 
compensatory proliferation, in a linearor'‘onc-hit" fash¬ 
ion. For example, this may describe the action of a single 
hormone or growth (actor molecule interacting with a 
highly specific cell-surface receptor molecule that func¬ 
tions as a highly sensitive signal transducer. Although it 
is difficult to imagine how this type of response might 
also apply to relatively nondisti net molecules like chloro- 
methanes or their reactive metabolites, the possibility is 
suggested hy the ability of CCU. CHCI,. and CH,CI ; to 
act as activators of protein kinase C in a cell-free in vitro 
system (Rogbani ct a!.. ! 987} and by the mitogenic activ¬ 
ity of even very small concentrations of CCU to human 
embryonic lung fibroblasts growing in culture (Brun'tus. 
1987). However, results obtained with cell-free systems 
and in vitro systems involving proliferating embryonic 
cells may have no relevance to the proliferative behavior 
of adult cells in vivo, one reason being that three-dimen¬ 
sional intercellular interaction is likely to play an impor¬ 
tant regulatory role in adult nonregenerating tissue. 

The chemical induction of cell death with definitively 
or approximately single-hit kinetics is known to occur, 
but it generally involves DNA damage resulting tn a le¬ 
thal mitotic event due to replication prior to repair of 
that damage (i.e., repair that is either total or sufficient 
to allow replication to proceed). This has been demon¬ 
strated for benzof a] pyrene in cultured mouse liver epi¬ 
thelial cells (Landolph et at.. 1976) and fora number of 
alkylating agents and DNA-crosslinking agents in bacte¬ 
ria and DNA repair-deficient CHO cells (Thompson et 
at.. 1983: Hoy et at.. 1985). The compounds CCU, 
CHCI). and CH,C1 2 . however, are only very weakly ge- 
notoxic. if at all. to mammalian ceils in vivo. Further¬ 
more. even (O-week exposures of rats to the potent alkyl¬ 


ating agent dtethylnitrosamine failed to induce 
(apparently necrosis related) hepatocellular proliferation 
in strict proportion to the administered dose; rather, in¬ 
duced replication was found to have a quite nonlinear 
dose response (Deal et at, 1989). Extremely efficiently 
induced cell death, suggestive of single-hit kinetics appar¬ 
ently not mediated by DNA damage but rather by dam¬ 
age to the plasma membrane, has been observed in chem¬ 
ical- and radiation-exposed immature primordial oo¬ 
cytes from rodents and (even more dramatically so) 
from squirrel monkeys, but this greatly enhanced sensi¬ 
tivity is absent in growing or mature oocytes (Dobson 
and Felton. 1983). To the extent that this type of ultra¬ 
sensitive msponse is known to involve only highly spe¬ 
cialized immature cells, it would appear not to be rele¬ 
vant to hepatocytes under a lifetime, chronic exposure 
scenario. 

There is currently no evidence that CCU, CHCU, or 
CH-Cf can exert a proliferative effect on hepatocytes (or 
any other cell type) at extremely low doses in vivo. On 
the contrary . CCU and CHCI, are agents classically used 
to induce experimental hepatocellular necrosis and sub¬ 
sequent regenerative hyperplasia (Leevy et a/.. 1959: 
Klassen and Plaa. 1966: Reynolds. 1967; Schultze el a/.. 
1973: Moore et at.. 1982). Perhaps the most is known 
regarding the mechanism of CCU-induced hepatocellu¬ 
lar toxicity. The events preceding CCI 4 -induced cell 
death are believed at present to involve the creation of a 
reactive trichioromethyl free-radical metabolite predom¬ 
inantly within the endoplasmic reticulum (ER) by the 
action of microsomal cytochrome P450 enzymes—bul 
also by P450 activation at other cellular sites and by non- 
enzymatic activation by the electron transport chain 
within mitochondria—followed by secondary reactions 
of this metabolite with proteins and phospholipids that 
then interfere with macromoieeule synthesis and func¬ 
tion within the cell (Cheeseman et al.. 1985; Tomasi et 
at., 1987; Villamiel er a!.. 1987). One critical step in this 
process is now thought to be a resulting rapid inhibition 
of the ER CA 2+ pump, leading to leakage of sequestered 
CA 2t from the ER and, if severe enough, to a postulated 
lethal disruption in intracellular calcium ion concentra¬ 
tions (Moore, 1980; Long and Moore. 1986; Ray and 
Moore, 1986). It is interesting in this regard that altered 
ER Ca 2t pump activity was not found to be linearly re¬ 
lated to cytosolic CA-* concentration and that the latter 
was highly nonlineariy related to applied CCU (and even 
more so for applied CHCL) concentrations in studies us¬ 
ing primary cultures of rat hepatocytes (Long and 
Moore, 1986). In particular, this study-showed that while 
ER CA 2+ pump activity, 5 min after exposure of cells to 
0.3 mM CCU, underwent a significant decrease to two- 
thirdsofthat of nonexposed cells, the cytosolic Ca J+ con¬ 
centration (as measured by phosphorylase a activity) was 
not significantly increased, indicating that a quite size¬ 
able perturbation in ER Ca 2i pump activity was neces- 
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sary to disturb homeostasis in cytosolic Caconcentra¬ 
tion in the cells studied. 

It is, then, apparently unlikely that a single event at a 
single target mediates the type of CCI 4 -induced cell death 
that is known to stimulate hepatocellular proliferation; 
rather, in this case it is likely that the conjunction of mul¬ 
tiple events at some point constitutes a level of irrevers¬ 
ible cell injury that inevitably leads to eventual cell death 
(Smucklcr and James. 1984). Although probably multi¬ 
ple, such events mignl nevertheless all impinge on a sin¬ 
gle most critical cell function affected I possibly. e.g„ the 
maintenance of Ca : * homeostasis by the combined ac¬ 
tivity of ER Ca 2+ pump enzymes). Whether single or 
multiple, the disfunctions that mediate CCL-induced 
ceil deaih are likely to involve, at some stage, the satura¬ 
tion of one or more critical enzymes. Enzyme saturation 
can explain threshold-like dose-response phenomena in¬ 
volving even the simplest types of chemical reaction 
(Licko. 1971), The assumption that a given response, 
such as death of an individual animal orofan individual 
cell, requires the joint occurrence or interaction of 
exceedances of many such thresholds, and that these 
thresholds have a normal (Gaussian I distribution with 
respect to dose, gives nse to the lognormal (Or "log- 
profcif) dose-response model, which is one of the most 
widely used models i.n classical toxicology and pharma¬ 
cology for the description and analysis of dichotomous 
quantal dose-response data (Finney. 1957; Levine. 
1983 ). 

For the reasons discussed above, the relationship be¬ 
tween estimated values of acute effective dose 0 and cor¬ 
responding predicted values of the percentage of hepato- 
eytes in S-phase i %S) was modeled urine the lognormal 
dose-response model. 



-r %So. 


(J) 


in which T is the cumulative normal distribution func¬ 
tion. %5o is the background level of%5( i.e., ar p - 0 j, 
u is the location parameter (the dose corresponding to 
the median or geomelric mean response), and i r is the 
shape parameter (the geometric standard deviation of 
the lognormaily modeled response). The factor oi l00%/ 
3 and the corresponding assumed upper bound of (33rib 
+ %So) lor %.S used in Eq. (1) account lor the fact thal 
H>S represents a fraction of the totai hepatocytes under¬ 
going replicative DNA synthesis at 24 or 48 hr postexpo¬ 
sure, whereas not all hepatocytes induced into 5-phase 
are involved in synthesis precisely at that time. Clearly, 
a synchronous induction of 100% of hepatocytes into S- 
pbase would be lethal. Detailed experiments have shown 
that CCI 4 intoxication or partial hepateciomy in rats and 
mice induces a wave of regenerative hepatocellular pro¬ 
liferation corresponding to an S-phasc pulse commenc¬ 
ing between 24 and 48 hr alter administration and per- 
sisting between approximately 24 and 36hr( Leevy tHat. 


1959; Gerhard. 1973. Rabes el ai. 1973; Schultze et al.. 
1973). Such proliferative response kinetics have been 
shown to be quite heterogeneous within rat liver, with 
hepatocytes in most liver zones reaching maximum %$ 
values of between 20 and 35% at various times after in¬ 
duced acute liver injury (Gerhard. 1973; Rabes et at, 
1973). Furthermore, the average duration ol'S-phase in 
ratsisubcutd hr(Rabes, 1973), which isaboutone-third 
of the average (24-hr) duration of dramatically increased 
S-phase activity after chemically or physically induced 
hepatocellular trauma (Leevv et al.. 1959; Gerhard, 
1973; Rabes etui.. 197.3: Schultze era/.. 1973). This sug¬ 
gests that I - or 2-day interruptions in a chronic oral dos¬ 
ing regimen, such as the 5-day/week oral dosing regimen 
typical in cancer bioassays, would not substantially 
change the time-weighted average value of%> .5 associated 
with that regimen in comparison with that ofan uninter¬ 
rupted regimen. A nonresponse in average %Sto such an 
interrupted dosing regimen was therefore assumed in the 
present analysis for the compounds (CCL. CHCJ : ,.I con¬ 
sidered here that were orally administered in fcioassays 
on a 5-cay.-' week schedule. 

The lognormal model described above was fit to the 
daia shown in Tabic 2 using least-squares optimization 
of values for ^ and c, yielding the best-fit values listed in 
that i able, i Confidence bounds tor the lognormal param¬ 
eter estimates listed in Table 2 arc omitted—see Discus¬ 
sion.j The til ra the model to data on induction of9b S in 
hepatocytes of mice orally dosed (either once, or twice 
separated by 24 hr) with CCl d is shown in Fig, ]. 

The Doolrtilc et at (19S7) study included some i nves- 
tigation ol how chronic (daily) dosing affects %.Vand se- 
rdr; cfliyu.e icveL u mice as compared to single dosing. 
Their results indicate that induced %S is lower with 
chronic as Opposed to single CC1 4 Hosing by a factor of 
2 to 3. For the present analysis, it was assumed for all 
com pounds considered that the factor l00%/3 in Eq. (1) 
must be replaced by l0O%/9 when this model is applied 
to chronic dosing situations. The latter assumption im¬ 
plies that no chronic etfeclive dose can produce a value 
of % .9 greater than , where 


100% 

9 


%S a . 


( 2 ) 


Using this model, the predicted value of%5' nttl equals 
11 .6% for the Doolittle et al. (1987) study. 


CKAS Models Applied to Liver Cancer Data for Mice 

A mathematical description of the approximations for 
the two-stage exponential and geomelric CKM models 
used in the present analysis is given by Bogen ( I9S9). 
These models predict the expected value of cumulative 
cancer risk, R , as a function of age I minus a latency pe¬ 
riod here assumed to be 0) for a constant, lifetime expo- 
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TABLE 3 

Fitted Values for Growth and Mutagenic Potency Parameters in CKM Models” 
of Liver Tumorigenicity of Chlorinated Methanes to Mice 


Compound 

Assumed 
turn on genic 
mechanism 



CKM model type* 




Exponential 


GwttKtric 

g\ * gj 
(year -1 ) 

Si 

(year 1 ) 

b 

(Mg/kg liver)' 1 

Qy = ai 

(year -1 ) 

<32 

(year 1 } 

b 

(pg/kg liver)' 1 

CC1 4 

Promotion 

1.91 

>6.54 

0 

0.939 

>12,7 

0 


Mutation 

1.91 

1.91 

>1140 

0.939 

0.939 

>1090 

CHCI 3 

Promotion 

1.91 

6.18 

0 

0.939 

10.6 

0 


Mutation 

1.91 

1.91 

8.19 

0.939 

0.939 

7.79 

CH-Ci, 

Promotion 

2.51 

3.80 

0 

1.63 

3.93 

0 


Mutation 

2.51 

2.51 

1.74 

1.63 

1.63 

1.62 


“ See text for explanation of CKM dose-response models and their estimated parameters. The numerical subscripts 
, i on the growth parameters. ,£f, and a,. i = 1, 2, 3. refer to corresponding exposure endpoints I where exposure to 

| effective dose d is presumed to occur only during time / for r, i < t 2 . The parameter values listed are exact fits to 

dose-response data corresponding to control and low-dose bioassay responses. R 0 and R, ; , respectively (see Table t). 


sure scenario, using the relations given by Fqs.( tf-4) and 
(G-3) in the Appendix. To account for exposure scenar¬ 
ios involving n distinct exposure periods (i.e.. in cases 
involving a sequence of nsets of different parameter val¬ 
ues governing mutation and/or cell proliferation rales), 
the latter equations must be replaced by Eqs. (E-7 ) and 
(G-6), respectively, from the Appendix. All of the bioas¬ 
says considered here involved preexposure, exposure, 
and postejcpnsurt! periods. Therefore, to model hepato¬ 
cellular caneer risk to bioassay mice. Eqs. (E-7 ) and (G- 
6) were used with n set equal to 3 and with x(the expo¬ 
nent in the geometric CKM model) in £q. (G-6) set 
equal to 3. The three distinct exposure periods ending at 
times / - 1, 2, 3. respectively, used in each bioassay 
considered are shown in Table I. 

To lit the CKM models, it was assumed that bioassay 
mice have I.5X 10 8 hepatocytes/g liver tissue (Peterson, 
1985) at risk throughout their lives, of which 25% are 
diploid cells at risk for malignant transformation 
(Danielson, 1988). Assuming liver is 5.5% of mouse 
body weight ( U.S. EPA, 1988). the total number. N, of 
cells at risk is equal to body weight times 2 X 10 6 target 
cells/g body wt. It was further assumed that for nondosed 
animals the first and the second interstage mutation 
rates, m, -in 3 andM, = AT, (see Appendix (, are allequal 
to 5 X 10" 6 v" 1 (Bogen. 1989). The value of Kq 
= (A , )(tH|)(,Vf l ) for each bioassav data set considered is 
listed in Tabic 1- 

Next, the ceil growth parameters of the exponential 
and geometric CKM models, g and a. respectively (sec 
Appendix), were tit I exactly) to the data painsRe- Ra 


given in Table I for each bioassay under the assumption 
that the corresponding compound is a “pure promoter" 
within the CKM framework described above. Specifi¬ 
cally, in this first analysis, each of the compounds- was 
assumed to act only to increase the proliferation rate of 
cells in “initiated” or premalignant (stage-1) clones thai 
are themselves assumed to be-reristonT i., ity, cytotoxicity 
of necrogemc dosesbut nevertheless stimulated to prolif¬ 
erate to the same extent as are surrounding norma) hepa- 
tocytes during regenerative growth following dosing, 
which is consistent with experimental data on the com¬ 
parative proliferative behavior of ceils in rodent hepato¬ 
cellular foci and surrounding normal cells (Dragani et 
at., 1986). Here, d is the PBPK estimate of the lowest 
effective bioassay dose (also listed in Table I) to liver 
tissue at which a significantly increased liver tumor inci¬ 
dence was observed, R 0 is background tumor risk, and 
R d is observed tumor risk at effective dose d. Under tile 
pure promoter assumption, the somatic mutation rates 
during dosing, m 2 and M z , are equal to those for non¬ 
dosed animals. For each compound, two unique values 
were obtained for the growth parameters g and a corre¬ 
sponding to nondosed periods (g { = g 3 and u, =■ a j, 
respectively) and the dosing period (g 2 and U }. 
respectively) of the bioassav serving as the basis for risk 
extrapolation. The values obtained are listed in Table 3. 
The corresponding exponential and geometric CKM 
models fit to the two (control and low-dose) bioassav 
data points forCCL, are shown in Fig. 2a. 

To proceed with the "promotion-only” analyses, it 
was assumed that an increase in % 5 over %5i, by a given 
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FlO. 2. Two-Stage ccli-kinetic multistage (CKM I 
models (see Table 3 and corresponding text) fit to con¬ 
trol and low-dose bioassay data (see Table 1) on CC].- 
induced hepatocellular carcinoma incidence in female 
B6C3F1 mice. Dotted lines represent hypothetical 
CKM-predicted probability of tumor occurrence in the 
absence of proliferation of initiated (premaJtgnant) cells, 
(a) Promotion-only CKM models; (b) mutation-only 
CKM models; exponential and geometric models an- 
p.'uxmiaieiy coincide for dosed animals. 


fraction /of its predicted (or measured) bioassay-associ- 
ated range. Induced by a given effective dose D, causes a 
directly proportional increase in the net growth rate of 
premalignant cells. Thai is, it was assumed that 


r— y ^ So _ Z n — C Q 

%S d ~%So -d ~ Z»' 


%S 0 « %S l> ^%S ll ^%S m 
0 < .r < I 


0 «/« !, 


where z is a CKM growth parameter eq ual to either c or 
a\ the D subscript refers to a value corresponding to the 
Specified effective dose O: and the r/andO subscripts refer 


to values pertaining to dosed and control bioassay ani¬ 
mals, respectively. The upper bound of 1 on s given 
above accounts for the fact that ft.?is an indication only 
of cell replication rate, not of cell death rate. An en¬ 
hanced replication rate in proliferating premalignant 
ceils that is also associated with an enhanced cell death 
rate would lead to values of r less than 1. For the present 
analysis, in the absence of information on dose-response 
for compound-induced changes in cell death rates in 
vivo, it was assumed that s « 1. 

Analyses were next performed under the contrasting 
assumption that each of the compounds considered is 
only mutagenic and has no effect on premalignant cell 
proliferation, in these analyses, the first and second mu¬ 
tation rates during dosing, rr? 2 and .V/,, were assumed to 
be equal linear functions of effective dose £>, given by 

r»3 = M 2 = m 0 ( 1 + bD), (4) 

in which Pio represents the spontaneous mutation rate 
(m, = mi - M, - V/ 3 ) and h represents the relative muta¬ 
genic potency of effective dose D. Under this assump¬ 
tion, it follows, that gi = =gj and u, - a. = for each 
compound, and the corresponding calculated values of b 
are given in Table 3. (For example, the corresponding 
exponential and geometric CKM models fit To the bioas- 
sav data for CCU are shown in Fig. 2b-) To proceed with 
the "'mutation-only” analyses, the fraction /wasdefined 
simply as the ratio D/d for 0 «. D H d. 

Finally, in all analyses, background risk R 0 was as¬ 
sumed to be independent such that increased risk was 
defined as A - (R - R u )/{ i - R 0 ). where f? 0 K =s ^ 
and the vf and Osubs'-riptc csai.n refer to valuta, pertaining 
to dosed and control bioassay animals, respectively. 

The fitted CKM models, along with the assumptions 
stated above, were used to numerically calculate /as a 
function of A for each of the models and compounds 
considered. In each case, the particular, unique value of 
/denoted/*, was calculated to correspond to a predicted 
increased risk of,4 ' 10“/ These values are listed in Ta¬ 
ble 4. Next, the relations between D and /for the promot¬ 
ion-only and mutation-only CKM models were used to 
calculate the virtually safe effective dose D* correspond¬ 
ing to/* defined above. For the mutation-only analyses, 
the relation is simply/)* = (/*|(d). For the promotion- 
only analyses, the corresponding relation follows from 
Eq. (!), appropriately modified to reflect a chronic dos¬ 
ing assumption, and from Eq.(2),which together imply 
that 

D - 

For convenience, the magnitude of %S a relative to 
defined in Eq. (2) by shall be described using the 
ratio p defined as 
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TABLE 4 


Virtual Safe Dose" ofCtilorinaTED Methanes Extrapolated from Bioassay Data 
on mouse Liver Tlmors 


Compound 

Assumed 
lumongemc 
mechanism a 

Exponential CKM model* 

Geometric CKM model* 

Linearized 
multistage 
. model 4 

VSD 

(«g/day) 

/* 

(X10-‘) 

D* 

(mg/kg 

liver 1 

VSD 

(mg/day) 

/* 

(xicr*) 

D* 

(mg/kg 

liver) 

VSD 

(mg/day) 

ecu 

P 

>18 

<0.49 

<35 

>5.0 

<0.45 

<3.0 

_ 


M 

>2.4 

<5.2 

<2.5 

>2.4 

<5.2 

<2.5 

0.95 




X 10" 1 

x icr 4 


X lCr ! 

X HT 4 


CHClj 

P 

20 

170 

950 

6.3 ' 

150 

730 

— 


M 

7.2 

0.0071 

0.011 

3'.2 

0.007! 

0-0 il 

0.30 

ch 2 ci 2 

P 

9-0 

153 

9500 

6.1 

145 

9 J00 

— 


M 

5.0 

0.0060 

2.7 

5.0 

0.0060 

2.7 

23 


° Here, the virtual safe dose(VSD) is defined as that daily dcse ingested by a 60-kg human over a lifetime associated 
wish a predicted cancer risk of I0" 6 . 

°See texl for explanation ofCKM dose-response models and their estimated parameters. 
c VSDs listed for CCli, CHClj.and CH.Ch were calculated usingihe corresponding 95% upper-bound carcinogenic 
potencies of 0.063. 0.20. and 0.0026 (tng/kg/day)' 1 . respectively, derived by the U.5. Environmental Protection 
Agency (U.S. EPA, 1984. 1985a.b) using the “linearized" multistage dose-response extrapolation model {described 
in those references) applied to data on hepatocellular carcinoma induction in female B6C3F1 mice from all dose 
groups in the bioassays considered in the present study (see Table 2). The potency derived for CC1 4 was based on 
pooled data for male and female mice; that for CH.Cl, was based on data on incidence of hepatocellular carcinoma 
or adenoma. Metabolism was considered in potency estimation oniy in the ease ol'CHClj. where94% metabolism of 
the administered bioassay doses were assumed (U.S. EPA. 1985a). 
d P. promotion only! M. mutation only. 



whereupon it follows from Eqs. (3). (5). and (6) that 

. (7) 

For the promotion-only analyses, then, D * was calcu¬ 
lated using Eq. (7) in which /* was substituted for /, 
where it was assumed that t = 1 (see discussion above) 
and where the values for p, e. and p are listed in Table 2. 

Extrapolation of Liver Cancer Kisk to Unmans 

It »vas assumed that the virtually safe effective dose ob¬ 
tained for female bioassay mice applies to humans as 
well, i.e., that human hepatoevtes are as equally suscepti¬ 
ble to compound-mediated enhanced proliferation as 
mouse hepatocytes. To complete the CKM-based risk 
analysis. PBPK models for reference 60-kg (females) hu¬ 
mans were used with parameter values obtained as dis¬ 
cussed above to estimate a human daily effective dose. 
ZX as a function of a chronic ingested dose (c.g., daiiv 


uptake from drinking water) for each of the compounds 
considered. Exposures were presumed to occur once per 
day. each modeled to result in a constant 6-min infusion 
of 100% of ingested dose into the liver compartment, and 
simulated daily exposures were repeated ( for 7 to 10 sim¬ 
ulated exposure days) until virtual dynamic equilibrium 
in daily effective dose was attained. 

RESULTS 

Best-fit parameter values obtained for the 
lognormal models fit to data on induced cell 
proliferation for each of the compounds are 
shown in Table 2. The promotion-only and 
mutation-only CKM models fit to mouse 
bioassay data on hepatocarcinogenicity for 
each of the comopunds are summarized in 
Table 3. Based on the information from Ta¬ 
bles 2 and 3, as explained under Methods, 
calculated values of a. virtually safe effective 
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Fir,. 3. Hepatotoxic effective dose to liver tissue of a 
60-kg human (reference female) chronically exposed to 
a daily single ingested dose of CC1 4 (medium-dashed 
curve). CHCIj (short-dashed curve!, or CH.Cl; (solid 
curve), as predicted by physiologically based pharmaco¬ 
kinetic models described in the text. Effective dose was 
assumed to be total 1-hr metabolized dose postexposure 
(in mg/ ke liver tissue) for CCh and CHCI,. whereas for 
CHsCU it was assumed to be total daily dose (in mg/kg 
liver tissue) metabolized in liver by the glutathione -V- 
transferase pathway. 


dose. D + . for the different compounds and 
CKM models are listed in Table 4. 

The relations between chronic applied 
dose and a PBPK-predicted human effective 
dose fo: each of the compounds are shown in 
Fig. 3. This information was used to obtain a 
human applied VSD corresponding to each 
value of D* listed in Table 4. These calcu¬ 
lated VSDs for the different CKM models are 
shown in Table t, where they are compared 
to VSDs for eac h compound calculated using 
the ’‘linearized” multistage model applied to 
the same mouse bioassay data on induced 
hepatocarcinogenesis that served as the basis 
for CKM-based risk: extrapolation. 

DISCUSSION 

Because the CKM-based approach to can¬ 
cer nsk assessment illustrated here for chlori¬ 
nated methanes allows for a threshold-like 
dose response for agents considered to be 
“pure promoters,” corresponding VSDs asso¬ 
ciated with exposure to such agents will typi¬ 


cally be much lower than those based on non¬ 
threshold risk-extrapolation methods that 
presume increased cancer risk for these com¬ 
pounds is due solely to increased mutation 
rates that are linear as a function of dose at 
low doses. Thus, in this analysis for chlori¬ 
nated methanes, the VSDs calculated using 
promotion-only CKM models are between 
three and five orders of magnitude higher 
than those calculated using mutation-only 
CKM models and are from three to more 
than seven orders of magnitude higher than 
those calculated using the linearized multi¬ 
stage model. Using the promotion-only 
CKM models, the geometric model yielded 
slightly but consistently lower (more conserv¬ 
ative) VSDs than the exponential one. as ex¬ 
pected {Bogen 1989). 

It must be emphasized that many uncer¬ 
tainties remain in defining the most appropri¬ 
ate way to implement CKM-based cancer 
risk extrapolation, even assuming the CKM 
approach is preferable. The most fundamen¬ 
tal source of uncertainty lies in the assump¬ 
tion made regarding the extent to which the 
observed carcinogenicity of a given com¬ 
pound (e.g.. in bioassay animals) is due to 
promotion (enhanced cel] proliferation 
rates) alone, which in some cases might rea¬ 
sonably have a threshold-like dose-response 
relationship, and how much is due to in¬ 
creased mutation rates, which in two-stage 
CKM models have a linear-quadratic dose- 
response relationship. Clearly, if a compound 
is assumed to have any degree of mutagenic¬ 
ity at all, this effect will ensure that the CKM- 
predicted increased risk as a function of 
effective dose will be linear at sufficiently low 
doses (Portier, 1987). In particular, it is eas¬ 
ily shown that if, within a CKM framework 
in which induced cell proliferation is as¬ 
sumed to have a threshold dose response, in¬ 
creased mutation actually accounts for P per¬ 
centage of an observed increased tumor inci¬ 
dence at high doses, then increased risk 
calculated using a conservative mutation- 
only CKM model extrapolated to very low 


doses will c 
approxima 
Some of 
the foregoi 
for chlorin 

(1) the 
chronic do 
istered dos 

(2) the 
(compoun 
interspecie 
toxic (cell- 
on hepatoe 

(3) the s 
the CKM . 
of CKM ti 
geometric. 

(4) the i 
changes in 
changes in 
minus deal 

(5) the 
available o 
duced by 
CHiCb an 
range large 
as a logne 
response; 

(6) the r 
induced by 
long-term 
experience 

(7) the s' 
the presen 
confidence 
rameter va 
the tumor 
groups froi 
elusion of 
groups; ant 

(8) the 
cancer inc 
doses less t 
NCI (197( 
that com pc 


PM3006451751 


Source: https://www.industrydocuments.ucsf.edu/docs/gkxj0001 











doses will overestimate true risk by a factor of Some of these uncertainties can be reduced 
approximately 100% / or eliminated by the acquisition of the appro- 

Some of the other sources of uncertainty in priate experimental data (Nos. 1,2,4-6, and 
the foregoing CKM-based risk extrapolation 8, above). In particular, there is a clear need 
for chlorinated methanes include: to gather much more detailed data on prolif¬ 

erative response to chronic dosing, especially 
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(1) the effects of corn oil vehicle and 
chronic dosing on ingestive uptake of admin¬ 
istered dose: 

(2) the selection of the most appropriate 
(compound-specific) PBPK measure(s) of 
mterspecies equitoxic effective dose For cyto¬ 
toxic (cell-proliferative) and genotoxic effects 
on hepatocytes: 

(3) the selection, assuming the validity of 
the CKM approach, of an appropriate form 
of CKJVt model to use, e.g„ exponential vs 
geometric, semistochastie vs fully stochastic: 

(4) the quantitative relationship between 
changes in measured increases in %.S' and 
changes in actual net proliferation (i.e.. birth 
minus death) rates in premalignant cells; 

(5) the very limited dose-response data 
available on increased in mouse liver in¬ 
duced by acute exposure, particularly for 
CHjCh and in general for doses spanning a 
range large enough to validate a model (such 
as a lognormal model) for low-dose dose 
response; 

(6) the relationship between increased %S 
induced by acute exposure to that induced by 
long-term chronic exposure similar to that 
experienced by bioassay animals; 

(7) the statistical uncertainty omitted from 
the present analysis by failing to consider 
confidence bounds on estimated CKM pa¬ 
rameter values and by using information on 
the tumor response of only the low-dose 
groups from the mouse bioassays to the ex¬ 
clusion of available data from higher dose 
groups; and 


among cells within putative preneoplastic 
foci in animals receiving very low chronic 
dose levels. Defining rather precise dose-re¬ 
sponse relations for such induced prolifera¬ 
tion should be experimentally achievable be¬ 
cause it is possible to score a sufficiently large 
number of cells (from even a relatively small 
number of experimental animals) to obtain 
statistically reliable dose-response character¬ 
izations (particularly with automated cell 
scoring techniques). Another uncertainty 
source (No. 7) is fairly easily addressed by the 
incorporation of standard statistical proce¬ 
dures within the general approach to CKM- 
based cancer risk extrapolation outlined 
above. This was not undertaken in the pres¬ 
ent investigation because its focus was on an¬ 
alytic and methodological issues concerning 
the distinction between the mutation-only 
and promotion-only approaches, as plausible 
extremes, to CKM-based cancer risk predic¬ 
tion. (However, it is unlikely that the wide di¬ 
vergence in risks predicted here using these 
two approaches would have been greatly re¬ 
duced had such statistical considerations 
been addressed.) Other sources of uncertainty 
(No. 3) are unlikely to be significantly re¬ 
duced until a better understanding of sponta¬ 
neous and chemically induced carcinogenesis 
is obtained. 

Despite the difficulties that are involved in 
addressing such uncertainties, the CKM ap¬ 
proach incorporates the emerging facts re¬ 
garding carcinogenesis to a far greater extent 
than do the more traditional cancer risk ex¬ 
trapolation models. Even when data avail- 



(8) the unknown dose response for liver able to implement a CKM-based risk extrap- 
cancer induction in mice receiving CCI 4 olation for a given compound are limited, 
doses less than the low-dose level used in the consideration of this quantitative approach 
NCI (J976a,b, 1977) cancer bioassays for may yield readily testable hypotheses. For ex- 
that compound. ample, if a suitable chronic dose level were 
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selected, the CKM initiator-only and pro¬ 
moter-only models for CHC1 3 -induced hepa- 
tocarcinogenesis in mice discussed above 
would predict very different tumor inci¬ 
dences in 1 O 0 exposed bioassy mice com¬ 
pared to those in controls, tn particular, if 30- 
g female B6C3F1 mice were administered a 
chronic daily dose of 67 mg/kg CHCI 3 by ga- 
vage, the PBPK model for this compound 
used earlier would predict a corresponding I- 
hr effective dose of 1000 mg/kg liver tissue. 
This effective dose, in turn, would induce a 
predicted increase in the value of the cell 
growth parameter, in either the exponential 
or geometric CKM protomoter-onlv models, 
equal to 10 % of that necessary to raise the in¬ 
cidence from its spontaneous rate of 1.25% to 
the rate ot 88.9% observed in bioassay mice 
given a chronic oral dose of 230 mg/kg/day 
(NCI, 1976c; Reuber, 1979; see Table 1), ac¬ 
cording to the lognormal model of induced 
cell proliferation for CHQ 3 specified in Table 
2. This proliferative response, in turn, would 
be predicted to result in an incidence rate of 
hepatocellular carcinoma between 1.9 and 
3.7%, according to the CKM promoter-only 
modek fot-CHC'-j specified in Tabic 2. How¬ 
ever, according to the corresponding muta¬ 
tion-only models, the incidence rate should 
be 41%. Time-to-tumor information from 
bioassays involving serial terminations 
would, of course, greatly improve the ability 
to test a promoter-only hypothesis within the 
CKM framework and-, based on spontaneous 
tumor incidence data alone, might even be 
used to test exponential CKM models against 
their geometric counterparts. 

Finally, it should be emphasized that even 
if the CKM framework accurately character¬ 
izes both spontaneous and chemically in¬ 
duced carcinogenesis, some—perhaps even 
most—compounds experimentally observed 
to be carcinogenic to animals are likely to 
cause this response by eliciting a mixture of 
mutagenic and cell-proliferative processes. 
There is some degree of correlation between 
mutagenicity and cytotoxicity (Grisham and 


Smith, 1984; Peterson and Peterson, 1985). 
It may therefore be that, as a consequence of 
associated cytotoxicity, some classically mu¬ 
tagenic carcinogens indirectly elicit a prolifer¬ 
ative response sufficient to account for some 
or most of their observed carcinogenicity 
(within the CKM framework), as suggested 
by the recent data of Deal et al. (1989) on 
diethylnitrosamine-induced hepatocellular 
proliferation in rats. It may also be that carci¬ 
nogenic compounds with little observed ge- 
noto.xic activity in mammalian cells in vivo, 
such as CC1 4 , may exert their carcinogenic 
effect solely by means of mutation (perhaps 
by some as yet unidentified mechanism). Be¬ 
cause mutagenic potential can never be ruled 
out definitively, the most controversial issue 
involved in the application of CKM models 
to cancer risk assessment is likely to be the 
definition of the necessary elements of any 
data base sufficient to warrant a reasonable 
conclusion that a given compound is very un¬ 
likely to be genotoxic to mammalian cells in 
vivo, such that it might be considered a pure 
promoter within the CKM framework. 

APPENDIX 

For CKM models in general, the stage-,/ 
cell populations of a A-stage process of carci¬ 
nogenesis are presented by the stochastic pro¬ 
cesses Xj{ T),j~ 1,.... k, over the period 0 

t < T, where 7t}(0) = 0: X 0 ( 7 ) is a deter¬ 
ministic nonnegative (DN) (“driving") func¬ 
tion; and for / - 0. k - 1 and 

m ; ( I ) T< 1, are DN functions (stage-specific 
mutation rates) such that the probability of a 
unit gain by X J+ \(T) caused by a simulta¬ 
neous unit loss from Xj(T) during the inter¬ 
val (T,T+ dT) is approximately m,{ T)dT , 
and the probability of more than one such 
event occurring in this interval is negligible. 
In semistochastic CKM models, uy( T) J - l, 

■ ■ ■, k — I, are DN impulse-response (or 
weighting) functions used to represent 
(premalignant) cell proliferation subsequent 
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to the occurrence of each new stage-; clone 
arising by mutation of a stage-(; - 1) cell, 
where w k (T) = 1. Finally, let R(t) 

= Prob { X k (t ) 3= 1}, the cumulative distribu¬ 
tion of the waiting time for the occurrence by 
time t of the first increment in X k (t). repre¬ 
senting the occurrence of the first cancer cell, 
which here is presumed to result rapidly in a 
lethal tumor. It has been shown (Bogen, 
1989) for semistochastic CKM hiodels in 
general that R(t) ~ 1 — exp(£A",t(r)), where 
E denotes the expectation operator. 

Below are the approximate solutions for 
two-stage (i.e., k = 2), semistochastic, ex¬ 
ponential, and geometric CKM models of 
carcinogenesis, for the case where parameter 
values governing mutation rates and cell pro¬ 
liferation in intermediate-stage (here, in 

stage-1) cells change at times t, for i = 1. 

n , where t 0 = 0, l,-.i « and t„ = t. Here it is 
assumed that 2f 0 ( T) is equal to the constant 
X„ and that m^lj- \ < T « t,) for j - 0 and 
j = I are equal to the constants in, and .V/,, 
respectively. From these assumptions, it fol¬ 
lows (Bogen, 1989) that for both exponential 
and geometric CKM models. 


x l {t,) = ,y 0 —(^' r, '- |) - 1) (E-l) 
g, 

XjU,)=X i o 

si 

-gAt.-ti- i)~l) (E-2) 

and note that, because Wj ( T) = e gT when n 
= 1, Eqs. (A-l) and (A-2) imply that 

(E-3) 

and 

(E-4) 

To consider cases in which n> 1, we now de¬ 
fine 

f ‘ 

X g;{t, - l,-i) forrKs 
O., = 1 (E-5) 


.0 


for r> s 


and note that vv,( T) = exp(g,T + dur¬ 
ing the period r„_, to t„ for stage-1 clones aris¬ 
ing at time T where t,-, t,. Upon substi¬ 

tution and integration, Eqs. (A-l) and (A-2) 
yield 


EX,(t„) 


Z J x t 

i-1 1 


and 

EXzUn)** EX 2 (t„^) 


f" 

J t*-t 


EX,(t„) « 2 A'i(ii)exp(f, + l ,J 1 (E-6) 
/=? 

and 

n 

EX ,( T)dT. EX 2 {t„) ** 2 X 2 (t,) 


X 0 m-! vv, ( t, — T)dT. 

(A-l) 


(A-2) 

In the following evaluations of Eqs. (A-l) 
and (A-2), the growth response function. 
W\(T), for all proliferating premalignant 
(i.e.. stage-1) cells is taken to be a function of 
the growth parameter z, where z is designated 
g and a in exponential and geometric CKM 
models, respectively, and where z at time T 
equals the constant z, for <£ T < t h i - !. 
... , n. 

/. Exponential Two-Stage CKM Models 

For two-stage exponential CKM models, 
we first define the functions 


+ 2 tj)Xz'expU t ^.,. t ). (E-l) 


[I. Geometric Two-Stage CKM Models 

For two-stage geometric CKM models in 
the case where n = 1, it has been shown 
(Bogen, 1989) that w,(7') = (aT + I) 1 
(where, e.g., jc = 3 if premalignant clones are 
assumed to grow spherically with radial 
growth ultimately proportional to time). For 
these geometric CKM models, we first define 
the function 
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m t M. 


X 1 (l,)=X 0 - 

°(x+ l)(j:+2)a? 

X - T-i) + 1 )' r+: 


-(Jf+2)a,(/,-r w )+ 1) (G-l) 

and note that from Eqs. (A-1) and (A-2) it 
follows that 


EXdt ' > ~ x » r~frr [(«■<■ + 1 ) x * ' - 1 h 

(-V+ 1 )<J| 


(G-2) 


and 


EX 2 ( h ) ~i 2 {r,): (G-3) 

To consider cases in which n > 1, we now de¬ 
fine 


I ! + Z u >(t, ~ h-\) for r < s 

7oH (G-4) 

II for r > s 


and note that w, ( T) = ( a. T -F y ,during 
the period r„_., to f„ for stage-1 clones arising 
at Time 7 where t,_, ^ T i,. Upon substitu¬ 
tion and integration, Eqs. (A-l) and (A-2) 
now yield 


EX i ((„) —. 2^ Ao - ■ — 

, {x + 1 )a t 

and 

t.x 2 {t„) sb 2 xt 2 Ui) 

I-1 


(G-5) 


+ 2 


x HiMl 

[ 0 (,x + 1 )(x + 2 )a t a, 


V 2 _ -. x+ Z _ _,x-f2 I x+2 x 

vT t.j Ti.j- 1 Ti'+I.y ' 7 /+) J- ) } 


(G-6) 
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